Three-Dimensional  Hydrodynamic  Model  for  Prediction  of 
Falling  Cylinder  Through  Water  Column 


i 


Peter  C.  Chu,  chu@nps.navy.mil 
Chenwu  Fan,  fan@nps.navy.mil 
Ashley  Evans,  ashley.evans@navy.mil 
Anthony  F.  Gilles,  anthony.gilles@navy.mil 
Naval  Postgraduate  School,  Monterey,  CA  93943 
Peter  Fleischer,  fleischerp@navo.navy.mil 
Naval  Oceanographic  Office,  Stennis  Space  Center,  MS  39529 


i 

Abstract — A  three  dimensional  hydrodynamic  model 
based  on  triple  coordinate  systems  is  developed  to  predict 
translation  and  orientation  of  falling  rigid  cylinder  through 
the  water  column:  earth-fixed  coordinate  (E-coordinate), 
cylinder’s  main-axis  following  coordinate  (M-coordinate), 
and  hydrodynamic  force  following  coordinate  (F- 
coordinate).  Use  of  the  triple  coordinate  systems  and  the 
transforms  among  them  leads  to  the  simplification  of  the 
dynamical  system.  The  body  and  buoyancy  forces  and  their 
moments  are  easily  calculated  using  the  E-coordinate 
system.  The  hydrodynamic  forces  (such  as  the  drag  and  lift 
forces)  and  their  moments  are  easily  computed  using  the  F- 
coordinate.  The  cylinder’s  moments  of  gyration  are  simply 
represented  using  the  M-coordinate.  A  recursive  model  is 
developed  on  the  base  of  the  triple  coordinate  transform  to 
predict  the  cylinder’s  translation  and  orientation.  To 
evaluate  the  recursive  model,  a  cylinder-drop  experiment 
was  conducted  at  the  Naval  Postgraduate  School  swimming 
pool  in  June  2001.  Comparison  between  the  modeled  and 
observed  results  shows  a  great  potential  of  using  the  triple 
coordinate  transform. 

Key  Words —  Triple  coordinate  systems,  momentum  of 
cylinder,  moment  of  momentum  of  cylinder,  cylinder  drop 
experiment  (CIDEX). 

1.  INTRODUCTION 

Study  on  the  movement  of  a  rigid  body  in  fluid  has  wide 
scientific  significance  and  technical  application.  The 
theory  of  dynamics  of  a  rigid  cylinder  allows  one  to  set 
up  six  nonlinear  equations  for  the  most  general  motion: 
three  momentum  equations  and  three  moment-of- 
momentum  equations.  The  scientific  studies  of  the 

1  This  work  was  sponsored  by  the  Office  of  Naval  Research 
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hydrodynamics  of  a  rigid  cylinder  in  fluid  involve  the 
nonlinear  dynamics,  flight  theory,  body-fluid  interaction, 
and  instability  theory.  The  body  forces  include  the 
gravity  and  the  buoyancy  force.  The  hydrodynamic 
forces  include  the  drag  and  lift  forces  that  depend  on  the 
fluid-to-body  velocity.  Usually,  a  nonlinear  dynamical 
system  is  needed  to  predict  a  falling  rigid  body  in  fluid, 
e.g.,  von  Mises  [1]. 

When  the  rigid  body  is  cylindrical,  this  dynamical 
system  can  be  simplified  using  three  coordinate  systems: 
earth-fixed  coordinate  (E-coordinate),  cylinder’s  main- 
axis  following  coordinate  (M-coordinate),  and 
hydrodynamic  force  following  coordinate  (F-coordinate). 
The  origin  of  both  M-  and  F-coordinates  is  at  the 
cylinder’s  COM.  The  body  forces  and  their  moments 
are  easily  calculated  using  the  E-coordinate  system.  The 
hydrodynamic  forces  and  their  moments  are  easily 
computed  using  the  F-coordinate.  The  cylinder’s 
moments  of  gyration  are  simply  represented  using  the  M- 
coordinate.  In  this  study,  a  recursive  model  is 
established  to  predict  the  cylinder’s  translation  velocity 
and  orientation  on  the  base  of  the  triple  coordinate 
transformation,  e.g.,  Crane  and  Duffy  [2]. 

To  evaluate  the  recursive  model,  a  cylinder-drop 
experiment  was  conducted  at  the  Naval  Postgraduate 
School  swimming  pool  in  June  2001.  It  consisted  of 
dropping  cylinders  into  the  water  and  recording  the 
position,  orientation,  and  time  using  two  digital  cameras 
at  (30  Hz)  as  the  cylinders  fell  2.4  meters  to  the  pool 
bottom. 

Theoretical  model  development  and  cylinder  drop 
experiment  are  depicted  in  this  paper.  The  outline  of  this 
paper  is  as  follows:  A  description  of  the  dynamics  of 
falling  cylinder  is  given  in  Section  2.  A  depiction  of  the 
model  cylinders  and  Cylinder  Drop  Experiment  is  given 
Sections  3  and  4.  Experimental  data  retrieval  and 
analysis  are  given  in  Section  5.  The  experimental  results 
are  described  in  Section  6.  Model-data  inter-comparison 
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is  presented  in  Section  7.  The  conclusions  are  listed  in 
Section  8. 

2.  DYNAMICS 


A.  Three  Coordinate  Systems 

Consider  an  axially  symmetric  cylinder  with  the  centers 
of  mass  (X)  and  volume  (B)  on  the  main  axis  (Fig.  1). 
Let  ( L ,  d,  X)  represent  the  cylinder’s  length,  diameter, 

and  the  distance  between  the  two  points  (X,  B).  The 
positive  X  -values  refer  to  nose-down  case,  i.e.,  the 

COM  is  lower  than  the  COV.  Three  coordinate  systems 
are  used  to  model  the  hydrodynamics  of  falling  cylinder 
through  the  water  column:  E-,  M-,  and  F-coordinate 
systems.  All  the  systems  are  three-dimensional, 
orthogonal,  and  right-handed.  The  benefit  is  to  simplify 
the  following  computations:  the  body  and  buoyancy 
forces  and  their  moments  in  the  E-coordinate  system,  the 
hydrodynamic  forces  (such  as  the  drag  and  lift  forces) 
and  their  moments  in  the  F-coordinate,  and  the  cylinder’s 
moments  of  gyration  in  the  M-coordinate. 


The  E-coordinate  system  is  fixed  to  the  Earth  surface 
with  horizontal  sides  as  x  and  y- axes  (along  the  two  sides 
of  the  pool),  and  vertical  direction  as  z-axis  (upward 
positive,  Fig.  2).  Suppose  a  cylinder  falling  into  the 
water  column.  The  cylinder  rotates  around  its  main  axis 
(ri)  with  an  angle  y/\  and  an  angular  velocity  of  £1 .  Its 
position  is  represented  by  the  center  of  mass  (COM),  and 
its  orientation  is  represented  by  two  angles:  y/2  and  y/3 

(Fig.  3).  Here,  y/2  is  the  angle  between  the  rraxis  and 
the  horizontal  plane;  and  y/3  is  the  angle  between  the 
projection  of  the  main  axis  in  the  (x,  y)  plane  and  the  x- 
axis.  The  angle,  y/2  +  7T  /  2,  is  usually  called  attitude. 

The  relative  coordinate  is  rigidly  connected  with  the 
cylinder.  The  origin  (O)  of  the  relative  coordinate  system 
coincides  with  the  center  of  mass  (COM);  the  rraxis  is 
along  the  central  line  of  the  cylinder;  the  r2-axis  is 
perpendicular  to  the  plane  constructed  by  rraxis  and  z- 
axis  (rrz  plane);  and  the  r3-axis  lies  in  the  (rrz)  plane 
and  is  perpendicular  to  rraxis.  The  selection  of  axes  (x, 
y,  z)  and  (r1?  r2,  r3)  follows  the  right-hand  rule.  Let  V*  = 
(Ci*,  C2*,  C3*)  be  the  three  components  of  the  velocity 
of  COM,  i.e.,  be  the  origin  velocity  of  the  coordinate 
system  (r1?  r2,  r3).  The  geometric  center  (GC)  is  located 
at  (j,  0,  0).  For  GC  below  COM,  x  >  0,  and  for  GC 
above  COM,  x  <  0.  The  relative  coordinate  system  (ru 
r2,  r3)  is  obtained  through  the  translation  and  two 
rotations  y/2  and  y/3  of  the  earth  coordinate  system.  Let 


the  position  vector  P  be  represented  by  PE  and  PB  in  the 
earth  and  relative  coordinate  systems,  PE  and  PB  are 
connected  by 
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where  (xm*,  ym*,  zm*)  represent  the  position  of  COM  in 
the  earth  coordinate  system. 

B.  Nonlinear  Dynamical  Equations 

Any  solid  object  falling  through  a  fluid  (air  and  water) 
should  obey  two  physical  principles:  (1)  momentum 
balance  and  (2)  moment  of  momentum  balance.  Let  Vw  = 

(Fwi,  Vw2,  Vw3)  be  the  water  velocity,  and  (co* ,co*,co*) 

be  the  components  of  the  angular  velocity,  referring  to 
the  direction  of  the  relative  coordinate  system.  The 
independent  and  dependent  variables  are  made  non- 
dimensional  by: 


where  g  is  the  gravitational  acceleration,  and  L  the  length 
of  the  cylinder.  The  non-dimensional  momentum 
equations  for  COM  are  given  by  (Mises  1959) 


dV  p-p  F 

— L  +  coV  -coV  = - -sin w  + — 1 — ,  (2a) 
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where  n  is  the  volume  of  the  cylinder,  p  the  cylinder 


density,  pw  the  water  density,  (iq*,  F2*,  F3*)  the 

components  of  water  drag.  The  non-dimensional 

equations  of  the  moment  of  momentum  for  axial 
symmetric  cylinder  are 
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where  x  *  is  die  distance  between  COM  and  geometric 
center  (GC),  (Mi*,  M2*,  M3*)  the  components  of  the 
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moment  due  to  water  drag,  (Ju  J2 ,  X3)  the  three 
moments  of  gyration, 

J,  =  {  (r22  +  r')dm*,  J2=  J  (n  +  r')dm*, 

+rl)dm*  .  (4) 

The  orientation  of  the  cylinder  ( y/2 ,  )  is  determined  by 

du /  dll r 

— L  =  o)2 ,  cos  i/s2 — ~  =  cor  (5) 

dt  dt 

The  eight  non-dimensional  nonlinear  equations  (2a-c), 
(3a-c),  and  (5)  are  the  basic  system  for  determining  the 
cylinder  movement  in  the  water  column.  Chu  et  al.  [3] 
developed  a  triple  coordinate  transform  method  to  solve 
the  dynamical  equations  (3)-(5).  The  model  outputs  are 
compared  to  the  experimental  data  collected  during  the 
Cylinder  Drop  Experiment  (CYDEX)  conducted  at  the 
NPS  swimming  pool  in  June  2002. 
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Fig.  1.  M-coordinate  with  the  COM  as  the  origin  X  and  (im, 
jm)  as  the  two  axes.  Here,  %  is  the  distance  between  the 

COV  (B)  and  COM,  (X,  d)  are  the  cylinder’s  length  and 
diameter. 


Figure  2.  Earth  coordinate  system. 
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Figure  3.  Cylinder  orientation  and  relative  coordinate 
system. 

3.  MODEL  CYLINDERS 

A.  Description 

Three  model  cylinders  were  used  for  the  drop 
experiment  at  the  Naval  Postgraduate  School  swimming 
pool.  All  had  a  circular  diameter  of  4  cm,  however  the 
lengths  were  15,  12  and  9  cm  respectively  (Chu  et  al. 
[5]).  The  bodies  were  constructed  of  rigid  plastic  with 
aluminum-capped  ends.  Inside  each  was  a  threaded  bolt, 
running  lengthwise  across  the  cylinder,  and  an  internal 
weight  (Fig.  4).  The  internal  cylindrical  weight  made  by 
copper  was  used  to  vary  the  cylinder’s  center  of  mass  and 
could  be  adjusted  fore  or  aft.  The  center  of  gravity  of  the 
model  cylinder  is  the  origin  of  the  body  fixed  coordinate 
system. 

The  model  cylinder  is  composed  of  six  uniform 
cylindrical  parts  (Fig.  4):  (a)  a  plastic  hollow  cylinder 
C(1)  with  mass  of  mi,  outer  and  inner  radii  of  Rx  and  R2 , 
length  of  (X  -2/!),  and  the  center  of  gravity  for  the  part 
(COMP)  to  be  at  its  geometric  center  located  along  the 
rraxis  at  r  =  x  ;  (b)  an  aluminum-capped  left  end  (Fig. 

5)  solid  cylinder  C(2)  with  mass  of  m2,  radius  of  Rx, 
length  of  /1?  and  COMP  located  along  the  rraxis  at  rx  = 
LI 2-  lxl 2  +%  \  (c)  an  aluminum-capped  right  end  solid 
cylinder  C(3)  with  mass  of  m3,  radius  of  Rx,  length  of  lu 
and  COMP  located  along  the  rraxis  at  rx  =  LI 2- lx/2  -  x ; 
(d)  a  cylindrical  thread  C(4)  with  mass  of  m4,  radius  of  X3, 
length  of  (L  -  21  x),  and  COMP  located  along  the  rraxis  at 
r,  =  x ;  (e)  a  cylindrical  threaded  bolt  C(5)  with  mass  of 
m5,  out  and  inner  radii  of  R2  and  R3,  length  of  l2 ,  and 
COMP  located  along  the  rraxis  at  rx  =  LI 2  -  x  -  h-  h!2\ 
(f)  an  adjustable  copper  cylindrical  weight  C(6)  with  mass 
of  m6,  outer  and  inner  radii  of  R2  and  R2,  length  of  /3, 
and  COMP  located  along  the  rraxis  at  rl  =  S  +  j, 
where  S  is  the  distance  between  the  COMP  of  the 
adjustable  weight  and  the  geometric  center  of  the  model 
cylinder. 
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B.  Moments  of  Gyration 

Since  the  six  parts  (all  cylinders)  all  have  uniform 
mass  distribution,  the  moments  of  gyration  for  these 
parts  are: 

Jl"  =  —(R*  +R2),  Jj  =  — - =—Rl 


Rl  -C  =y K ,  P  =  f-(R22+Rl), 

J'2"  =  =  —  [3£2  +  3/?2  +  (I  -  2/  )2  ] , 

3  12  1 

J(2)  =  J(2)  =  —(3R2  +  /2),  J(3)  =  J(3)  =  —  (3Z2  +  /2), 

12  1  1  2  3  12  ‘  1 

/^j-^[3^+(Z-2/)2], 
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J(5)  =  J(5)  =  —  (3Z2  +  3Z2  +/2), 

2  3  12  2  3  2 


J  ’  =JK  )  =  — (3 R  +3R~  +/  ), 

3  2  3  3* 


(6) 


where  the  superscripts  for  the  moments  indicate  the 
cylindrical  parts.  The  resultant  moments  of  gyration  is 
computed  by 


A  ,,,  2  L  — l  2  L-l 

J2=J,=ZjJ2  +m,z  +mS - -X)  +«( - +  J) 
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According  to  the  definition  of  COM,  the  coordinate  of 
GC  ( / )  is  determined  by 


[m  (L/2-/  -/2/2)-w/] 

z- 


(8) 


which  indicates  how  the  adjustable  weight  determines 
the  location  of  COM  for  the  model  cylinder. 


Figure  4.  Internal  components  of  the  model  cylinder. 
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Figure  5.  Internal  structure  of  the  model  cylinder. 


A.  Model  Parameters 
C.l.  Aspect  and  Density  Ratios 

Our  goal  was  to  choose  a  scale  that  was  somewhat 
representative  of  the  real  world  ratio  of  water  depth  to 
cylinder  length,  but  at  the  same  time  would  be  large 
enough  to  film  and  would  not  damage  the  pool's  bottom. 
The  model  cylinders  were  based  on  the  realistic 
assumption  that  a  3  m  cylinder  is  laid  in  water  depths  of 
45  m,  thus  producing  a  15:1  ratio.  The  depth  of  the  pool 
is  2.4  m.  From  this  ratio,  the  length  (Z)  of  the  model 
cylinder  is  chosen  as  15  cm.  The  addition  of  a  12  and  9 
cm  length  allowed  for  later  comparison  of  the  sensitivity 
of  water  phase  trajectory  to  the  ratio  of  cylinder  length 
over  diameter.  The  outer  radius  of  the  model  cylinder  is 
2  cm.  The  length/diameter  ratios  (L/D)  are  15/4,  12/4, 
and  9/4.  The  corresponding  density  ratios  (p/p  )  are 
1.70,  1.68,  and  1.88,  respectively. 


C2.  x  -Value 

In  each  of  the  three  model  cylinders,  the  location  of 
the  weight  (i,e.,  the  value  of  8 )  is  adjustable.  Use  of  (7) 
location  of  the  COM  (j  -value)  can  be  determined 
(Table  1).  The  positive  x  -value  indicates  that  COM  is 
below  the  geometric  center,  and  the  negative  X -value 
indicates  that  COM  is  above  the  geometric  center. 


Table  1.  Model  L/D,  density  ratio,  and  X -value  (unit:  cm). 


Cylinde 

r 

L/D 

pip. 

Zi 

2/2 

Zo 

X-\ 

X-2 

1 

15/4 

1.70 

1.85 

3.69 

0 

1.85 

3.69 

2 

12/4 

1.68 

1.21 

2.43 

0 

1.21 

2.43 

3 

9/4 

1.88 

0.68 

1.37 

0 

0.68 

1.37 

4 


5 


4.  CYLINDER  DROP  EXPERIMENT 

A  cylinder  drop  experiment  (MIDEX)  was  conducted 
at  the  NPS  swim  pool  in  June  2001.  The  purpose  of  the 
experiment  is  to  collect  data  about  cylinder’s  motion  in 
the  water  column  for  various  combinations  of  the  model 
cylinder  parameters.  It  basically  consisted  of  dropping 
each  of  three  model  cylinders  into  the  water  where  each 
drop  was  recorded  underwater  from  two  viewpoints. 
Figure  6  depicts  the  overall  setup.  The  controlled 
parameters  for  each  drop  were:  LI D  ratio,  x -value, 
initial  velocity  (Vin),  and  drop  angle.  The  Earth’s 
coordinate  system  is  chosen  with  the  origin  at  the  comer 
of  the  swimming  pool  with  the  two  sides  as  x-  and  y-axes 
and  the  vertical  z-axis.  The  initial  injection  of  cylinder 
was  in  the  (y,  z)  plane  (Fig.  2). 

A.  Initial  Velocity 

Initial  velocity  (Vin)  was  calculated  by  using  the 
voltage  return  of  an  infrared  photo  detector  located  at  the 
base  of  the  cylinder  injector.  The  infrared  sensor 
produced  a  square  wave  pulse  when  no  light  was 
detected  due  to  blockage  caused  by  the  cylinder’s 
passage.  The  length  of  the  square  wave  pulse  was 
converted  into  time  by  using  a  universal  counter. 
Dividing  the  cylinder’s  length  by  the  universal  counter’s 
time  yielded  Vin.  The  cylinders  were  dropped  from 
several  positions  within  the  injector  mechanism  in  order 
to  produce  a  range  of  Vin.  The  method  used  to  determine 
Vin  required  that  the  infrared  light  sensor  be  located 
above  the  water’s  surface.  This  distance  was  held  fixed 
throughout  the  experiment  at  1 0  cm. 

B.  Drop  Angle 

The  drop  angle  (initial  y/('n) )  was  controlled  using 

the  drop  angle  device.  Five  screw  positions  marked  the 
15°,  30°,  45°,  60°,  and  75°.  The  drop  angles  were 

determined  from  the  lay  of  the  pool  walkway,  which  was 
assumed  to  be  parallel  to  the  water’s  surface.  A  range  of 
drop  angles  was  chosen  to  represent  the  various  entry 
angles  that  air  and  surface  laid  cylinders  exhibit.  This 
range  produced  velocities  whose  horizontal  and  vertical 
components  varied  in  magnitude.  This  allowed  for 
comparison  of  cylinder  trajectory  sensitivity  with  the 
varying  velocity  components. 


Figure  6.  CIDEX  equipment. 


C.  Methodology 

For  each  mn  the  cylinders  were  set  to  a  x  -value. 
For  positive  x -value,  the  cylinders  were  placed  into  the 
injector  so  that  the  COM  was  located  below  the 
geometric  center.  For  negative  X  -value,  the  COM  was 
located  above  the  geometric  center  to  release.  A  series  of 
drops  were  then  conducted  in  order  of  decreasing 
cylinder  length  for  each  angle.  Table  2  indicates  number 
of  drops  conducted  for  different  drop  angles  and  X~ 
value  for  LID  =  15/4.  Number  of  drops  for  other  LID 
ratios  (12/4,  9/4)  is  comparable  to  that  for  LID  ratio  of 
15/4.  All  together  there  were  712  drops.  Each  video 
camera  had  a  film  time  of  approximately  one  hour.  At 
the  end  of  the  day,  the  tapes  were  replayed  in  order  to 
determine  clarity  and  optimum  camera  position. 

Table  2.  Number  of  drops  conducted  for  different  drop 
angles  and  X  -values  for  L/D  =  15/4. 


yV 

15° 

30° 

45° 

o 

O 

'sD 

75° 

13 

15 

15 

15 

12 

X\ 

9 

15 

15 

15 

9 

Z0 

12 

14 

15 

18 

6 

Z- 1 

0 

6 

6 

6 

0 

Z-2 

2 

6 

6 

0 

0 

5.  DATA  RETRIEVAL  AND  ANALYSIS 

A.  Data  Retrieval 

Upon  completion  of  the  drop  phase,  the  video  from 
each  camera  was  converted  to  digital  format.  The  digital 
video  for  each  view  was  then  analyzed  frame  by  frame 
(30  Hz)  in  order  to  determine  the  cylinder's  position  in 
the  (x,  z)  and  (y,  z)  planes.  The  cylinder's  top  and  bottom 
positions  were  input  into  a  MATLAB  generated  grid, 
similar  to  the  ones  within  the  pool.  The  first  point  to 
impact  the  water  was  always  plotted  first.  This  facilitated 
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tracking  of  the  initial  entry  point  throughout  the  water 
column.  The  cameras  were  not  time  synchronized;  thus, 
the  first  recorded  position  corresponded  to  when  the  full 
length  of  the  cylinder  was  in  view. 

B.  Source  of  Errors 

There  were  several  sources  of  error  that  hindered  the 
determination  of  the  cylinder’s  exact  position  within  the 
water  column.  Locations  above  or  below  the  camera’s 
focal  point  were  subjected  to  parallax  distortion.  Placing 
the  cameras  as  far  back  as  possible,  while  still  being  able 
to  resolve  the  individual  grid  squares,  minimized  this 
error.  Second,  the  background  grids  were  located  behind 
the  cylinder’s  trajectory  plane.  This  resulted  in  the 
cylinder  appearing  larger  than  normal.  This  error  was 
minimized  by  not  allowing  the  plotted  points  to  exceed 
the  particular  cylinder’s  length.  Third,  an  object  injected 
into  the  water  will  generate  an  air  cavity.  This  air  cavity 
can  greatly  affect  the  initial  motion,  particularly  at  very 
high  speeds  (hydro  ballistics).  The  air  cavity  effect  was 
deemed  to  be  minimal  due  to  the  low  inject  velocities 
used. 

C.  Data  Analysis 

The  3-D  data  provided  by  each  camera  was  first  used 
to  produce  raw  two  2-D  plots  of  the  cylinder’s  trajectory. 
Next,  2-D  data  from  both  cameras  was  then  fused  to 
produce  a  3-D  history.  This  3-D  history  was  then  made 
non-dimensional  in  order  to  generalize  the  results.  The 
non-dimensional  data  was  used  to  generate  impact  scatter 
plots  and  was  also  used  in  multiple  linear  regression 
calculations. 

6.  EXPERIMENTAL  RESULTS 

A.  Trajectory  Patterns 

After  analyzing  the  3-D  data  set,  seven  trajectory 
patterns  were  found.  The  plots  on  the  (y,  z)  plane  were 
chosen  for  trajectory  analysis,  as  this  plane  was  parallel 
to  the  direction  of  the  cylinder  drop.  The  generalized 
trajectory  patterns  are  described  in  Table  3,  and  Figures 
7-10.  The  water  phase  trajectory  a  cylinder  experiences 
ultimately  determines  the  impact  orientation.  In  MIDEX, 
the  categorizing  of  trajectories  into  general  patterns 
served  two  purposes.  Observed  trajectories  were  found 
to  be  most  sensitive  to  X  -value,  drop  angle  and  L/D 
ratio.  As  COM  distance  (z  -value)  increased  from  GC 
the  cylinder  tended  to  follow  a  straight  pattern.  As  COM 
was  moved  closer  to  the  GC  (decreasing  z -value)  the 
cylinder’s  trajectory  tended  towards  being  more  parallel 
with  the  pool’s  bottom.  At  steep  drop  angles,  the  cylinder 
experienced  little  lateral  movement  and  tended  towards  a 
straight  pattern.  Additionally,  as  LID  ratio  decreased 
more  complex  trajectory  patterns  developed.  This 
included  significant  oscillation  about  the  vertical  axis 
and  increased  lateral  movement. 

Table  3.  Description  of  relative  coordinate  based  trajectory 
patterns. 


Trajectory  Pattern 

Description 

Straight  (or  Slant) 

Little  angular  change, 

Nearly  parallel  to  z-axis, 

y/i  near  (90°  ±  15°) 

Spiral 

Oscillating  (no  rotation) 

Flip 

Only  once  as  z  *  <0 

Flat 

1/ /  near  0°  (no  oscillation) 

Seesaw 

x//^  near  0°  (oscillation) 

Combination 

Several  of  the  above  patterns 

Straight  Spiral  Flat 

L“  15cm;  Vi=  1  835m/s  L*  9cm:  1  553m/s  L=  15cm.  Vi=  218nVs 


Figure  7.  Trajectory  patterns. 

Flip  Seesaw  Flip-Flat 


L=  15cm:  Vi=2  874m/s  L=  15cm:  Vi=  1  379m/s  L«  9cm:  Vi»  2  294m/s 


Y<m>  Y  (m)  Y(m) 

Figure  8.  Trajectory  patterns. 
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COM  Position:  2 

Mine  Lenath: 

15 

12 

9 

Drop  Anqle:  15° 

Straiqht 

Straiqht 

Straiqht-Slant 

Slant-Straiqht 

Spiral 

Spiral 

Slant-Straiqht 

Slant-Straiqht 

Slant-Straiqht 

Slant-Straiqht 

Slant-Straiqht 

Spiral 

Drop  Anqle:  30° 

Slant-Straiqht 

Slant 

Spiral 

Straiqht 

Spiral 

Spiral 

Slant-Straiqht 

Straiqht 

Spiral 

Slant-Straiqht 

Slant-Straiqht 

Spiral 

Slant-Straiqht 

Slant-Straiqht 

Spiral 

Drop  Anqle:  45° 

Slant-Straiqht 

Spiral 

Spiral 

Slant 

Spiral 

Spiral 

Straiqht-Spiral 

Straiqht-Spiral 

Spiral 

Straiqht 

Straiqht 

Spiral 

Slant 

Slant-Straiqht 

Slant-Spiral 

Drop  Anqle:  60° 

Straiqht 

Straiqht-Slant 

Spiral 

Straiqht 

Straiqht 

Spiral 

Straiqht 

Straiqht 

Spiral 

Straiqht 

Straiqht-Spiral 

Straiqht-Spiral 

Straiqht 

Straiqht 

Spiral 

Drop  Anqle:  75  ° 

Straiqht 

Straiqht 

Spiral 

Straiqht 

Straiqht-Spiral 

Slant 

Straiqht 

Straiqht 

Spiral 

Straiqht 

Straiqht-Spiral 

Straiqht-Spiral 

Straiqht 

Straiqht-Spiral 

Straiqht-Spiral 

Figure  9.  Trajectory  patterns  for  %  =  %2. 


COM  Position:  -2 

Mine  Lenqth: 

15 

12 

9 

Drop  Anqle:  30° 

Flip-Straiqht 

Flip-Slant 

Flip-Straiqht-Spiral 

Flip-Straiqht 

Flip-Straiqht 

Flip-Straiqht-Spiral-Flip 

Drop  Anqle:  45° 

Flip-Straiqht 

Flip-Straiqht-Spiral 

Flip-Straiqht 

Flip-Straiqht 

Flip-Straiqht 

Flip-Straiqht-Spiral 

COM  Position:  -1 

Drop  Anqle:  30° 

Flip-Straiqht 

Flip-Slant 

Straiqht-Flip-Seesaw 

Flip-Slant 

Flip-Straiqht 

Flip-Straiqht-Spiral 

Drop  Anqle:  45° 

Flip-Spiral-Slant 

Flip-Slant 

Flip-Spiral 

Flip-Slant 

Flip-Slant 

Flip-Spiral 

Drop  Anqle:  60° 

Flip-Straiqht 

Flip-Straiqht 

Flip-Spiral-Seesaw 

Straiqht-Flip 

Slant-Flip-Slant 

Flip-Spiral 

Figure  10.  Trajectory  patterns  for  negative  COM 
position  ( x_\  ?  X-2  )•  Here,  COM  is  above  GC  when  the 
cylinder  enters  the  surface). 


B.  Impact  Attitude 

The  angle,  y/2  +  n  /  2 ,  is  the  impact  attitude  at  the 

bottom  of  the  water  column.  The  cylinder  burial  is 
largely  determined  from  the  impact  attitude  of  the 
cylinder.  Cylinders  whose  impact  attitudes  are 
perpendicular  (y/2=  90°)  to  the  sediment  interface  will 

experience  the  largest  degree  of  impact  burial  (Taber 
1999).  It  is  therefore  important  to  analyze  the 
relationship  between  impact  attitude  and  the  controlled 
parameters,  drop  angle,  Ein,  LID ,  and  X  •  The  experiment 
shows  that  both  LID  and  Kin  had  little  influence  on 
impact  attitude.  The  drop  angle  and  J,  however,  were 
the  largest  determinants  of  impact  attitude. 

From  Figure  11  it  is  apparent  that  there  are  several 
peaks  centered  near  90°,  140°,  and  180°.  Further  analysis 
reveals  that  these  peaks  correspond  to  the  COM  positions 

0,  1  and  2  (corresponding  to  j0  ,  X\  ,  X 2  )■>  respectively. 

COM  positions  -1  and  -2  (corresponding  to  X-\ ,  X-j) 

followed  the  same  trend  as  their  positive  counterparts. 

Although  drop  angle  was  not  the  most  influential 
parameter,  variations  did  induce  changes  in  impact 
orientation.  As  drop  angle  increased,  the  likelihood  of 
any  lateral  movement  decreased.  This  allowed  for  impact 
angles  that  were  more  vertically  orientated.  This  is 
primarily  due  to  the  fact  that  the  vertical  components  of 


velocity  were  greater  than  those  at  shallow  angles.  Thus, 
the  time  to  bottom  and  time  for  trajectory  alteration  was 
less. 


Figure  11.  Relationship  between  COM  position  and 
impact  attitude. 

7.  MODEL-EXPERIMENT  COMPARISON 

The  3D  model  predicts  the  motion  of  cylinder  inside 
the  water  column.  The  drop  angle  and  velocity  in 
CYDEX  are  not  the  initial  conditions  for  the  recursive 
model.  For  each  drop  in  CYDEX,  the  undersea  cameras 
measure  the  initial  translation,  velocity,  orientation,  and 
angular  velocity.  The  model-experiment  comparison 
starts  just  after  the  cylinders  are  fully  immersed  into  the 
water  column.  The  model  well  predicts  the  cylinder’s 
motion.  Here,  we  show  some  examples.  We  plan  to 
verify  the  model  using  NRL  data  (Holland  et  al.  [6]). 

7.1.  Positive  X  (Nose-Down) 

Example-1 

Cylinder  #1  ( L  =  15.20  cm,  ~p  -  1.69  gem’3)  with 
X=  0.74  m  is  injected  to  the  water  with  the  drop  angle 

45°.  The  physical  parameters  of  this  cylinder  are  given 
by 

J,  =  330.5  g  cm2,  J2=J3  =  5783.0  g  cm2 .  (9a) 

Undersea  cameras  measure  the  initial  conditions 
X0  =°’>;0  =°’Z0  =0’ 

uQ  =0,  v0  =  -1 .55  m  s’1 ,  wQ  =  -2.52  m  s’1 , 

^10  =  °>^20  =  60°>^30  =  “95° 

co10  =  0,  ©20  =  0.49  s"1,  ©30  =  0.29  s"1 .  (9b) 

Substitution  of  the  model  parameters  (9a)  and  the  initial 
conditions  (9b)  into  the  three-dimensional  hydrodynamic 
model  leads  to  the  prediction  of  the  cylinder’s  translation 
and  orientation  that  are  compared  with  the  data  collected 
during  CYDEX  at  time  steps  (Fig.  12).  Both  model 
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simulated  and  observed  tracks  show  a  slant- straight 
pattern. 

Example-2: 

Cylinder  #2  ( L  =  12.10  cm,  p  =1.67  gem’3)  with 
X  =  0.53  m  is  injected  to  the  water  with  the  drop  angle 
45°.  The  physical  parameters  of  this  cylinder  are  given 
by 

J]  =  271.3  g  cm2,  J2  =  J3  =  3206.5  g  cm2  .  (10a) 

Undersea  cameras  measure  the  initial  conditions 

xo  =  °’To  =  0,z0  =  0, 

uQ  =0,  v()  =  -1 .07  m  s’1 ,  w0  =  -1.74  m  s’1 , 

^10  =  °’^20  =  64°’^30  = -98°  ’ 

cow  =  0,  ffl20  =  4.50  s’1,  com  =  -0.21  s'1 .  (10b) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  CYDEX  at  time 
steps  (Fig.  13).  The  simulated  and  observed  tracks  show 
straight-spiral  pattern. 


Experiment  Model 


Fig.  12.  Movement  of  Cylinder  #1  ( L  =  15.20  cm, 
~p  -  1.69  gem3)  with  X =  0*24  m  and  drop  angle  45° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 


Experiment  Model 


Fig.  13.  Movement  of  Cylinder  #2  ( L  =  12.10  cm, 

_3 

p- 1.67  gem  )  with  x~  0.53  and  drop  angle  45° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 

Example-3: 

Cylinder  #3  ( L  -  9.12  cm,  p  =  1.88  g  cm’3 )  with  x  = 

0.58  m  is  injected  to  the  water  with  the  drop  angle  45°. 
The  physical  parameters  of  this  cylinder  are  given  by 

Jx-  235.0  gem2,  J2=J3  =  1556.8  gem2.  (11a) 

Undersea  camera  measure  the  initial  conditions 

xo  =  °>  T0  =  °>zo  =  °> 

uQ  =  0,  v0  =  -0.70  m  s’1 ,  wQ  =  -1.31  m  s’1 , 

W 10  =°^20  =  64°^30  =-90° 

O)l0  =  0,  co20  =  3.96  s’1,  O)30  =  0 .  (lib) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  CYDEX  at  time 
steps  (Fig.  14).  Both  model  simulated  and  observed 
tracks  show  a  spiral  pattern. 
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Experiment  Model 


Fig.  14.  Movement  of  Cylinder  #3  ( L  =  9.12  cm, 

_  _3 

/?  =1.88  gem  )  with  x~  0.58  m  and  drop  angle  45° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 

7.2.  Negative  X  (Nose-Up) 

Example-4: 

Cylinder  #1  (L  =  15.2  cm,  p  =  1.69  g  cm’3 )  with  x  = 
-1.48  cm  is  injected  to  the  water  with  the  drop  angle  30°. 
The  physical  parameters  of  this  cylinder  are  given  by 

Jx  -  330.5  g  cm2,  J2  =  J3  =  6233.8  g  cm2  .  (12a) 

Undersea  cameras  measure  the  initial  conditions 

xo  =  °’To  =  °>zo  =  °> 

uQ  =  0,  v0  =  -0.91  m  s’1,  w0  =  -0.78  m  s’1, 

^io  =  °>^20  =  17  V30  =  “96° 

colQ  =  0,  co20  =  -4.56  s’1,  co30  =  0.49  s’1 .  (12b) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  CYDEX  at  time 
steps  (Fig.  15).  The  simulated  and  observed  tracks  show 
flip-spiral  pattern.  The  flip  occurs  0.10  s  (0.12  s)  after 
cylinder  entering  the  water  in  the  experiment  (model). 
After  the  flip,  the  cylinder  spirals  down  to  the  bottom. 


Experiment  Madid 


_  -3 

p  =  1.69  gem  )  with  X =  -1.48  cm  and  drop  angle  30° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 

Example-5: 

Cylinder  #2  ( L  =  12.10  cm,  p  =  1.67  gem’3)  with 
X  =  -1 .00  cm  is  injected  to  the  water  with  the  drop  angle 
30°.  The  physical  parameters  of  this  cylinder  are  given 
by 

Jx  -  271.3  g  cm2,  J2=J3=  3312.6  gem2.  (13a) 
Undersea  cameras  measure  the  initial  conditions 
xo  =  °’To  =  °>zo  =  °> 

UQ  =0,  v0  =  -0.75  m  s’1,  WQ  —  —0.67  m  s’1, 

^10  =  °^20  =  24°’^30  =  "96° 

a>w  =  0,  co20  =  -5.08  s’1,  co30  =  0.15  s’1 .  (13b) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  CYDEX  at  time 
steps  (Fig.  16).  The  simulated  and  observed  tracks  show 
flip-spiral  pattern.  The  flip  occurs  0.11  s  (0.13  s)  after 
cylinder  entering  the  water  in  the  experiment  (model). 
After  the  flip,  the  cylinder  spirals  down  to  the  bottom. 
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EsiperimcfTl  MmJeE 


Fig.  16.  Movement  of  Cylinder  #2  ( L  =  12.10  cm, 

_  _3 

p  -  1.67  gem  )  with  %=  -1.00  cm  and  drop  angle  30° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 

7.3.  Zero  X  (Coincidence  of  COM  and  COV) 
Example-6: 

Cylinder  #2  (L  -  12.1  cm,  p  =  1.67  g  cm’3 )  with  x  = 

0  cm  is  injected  to  the  water  with  the  drop  angle  45°. 
The  physical  parameters  of  this  cylinder  are  given  by 

Jx  =  271.3  g  cm2,  J2  =  J3  =  3424.6  gem2.  (14a) 

Undersea  cameras  measure  the  initial  conditions 

xo  =  °’T0  =  °>zo  =  °> 

uQ  =  0,  v0  =  -1 .07  m  s’1 ,  w0  =  -1 .74  m  s’1 , 

^io  =  °^20  =  70V3o  =  -100° 

colQ>  =  0,  co20  =  -0.95  s’1,  co30  =  0.57  s’1 .  (14b) 

The  predicted  cylinder’s  translation  and  orientation  are 
compared  with  the  data  collected  during  CYDEX  at  time 
steps  (Fig.  17).  The  simulated  and  observed  tracks  show 
straight-seasaw  pattern.  However,  the  seasaw  pattern 
occurs  much  earlier  in  the  model  track.  Comparison 
among  Figs.  9-14  leads  to  the  fact  that  the  performance 
of  the  recursive  model  for  x  =  0  is  not  as  good  as  for 

Z*0. 


Experiment  MdcIcI 


Fig.  17.  Movement  of  Cylinder  #2  ( L  =  12.1  cm, 

_3 

p=  1.67  gem  )  with  x~  ®  cm  and  drop  angle  45° 
obtained  from  (a)  experiment,  and  (b)  recursive  model. 

8.  CONCLUSIONS 

(1)  Triple  coordinate  systems  are  suggested  to  predict  the 
translation  and  orientation  of  falling  rigid  cylinder 
through  water  column:  earth-fixed  coordinate  (E- 
coordinate),  cylinder’s  main-axis  following  coordinate 
(M-coordinate),  and  hydrodynamic  force  following 
coordinate  (F-coordinate).  The  benefit  is  to  simplify  the 
following  computations:  the  body  and  buoyancy  forces 
and  their  moments  in  the  E-coordinate  system,  the 
hydrodynamic  forces  (such  as  the  drag  and  lift  forces) 
and  their  moments  in  the  F-coordinate,  and  the  cylinder’s 
moments  of  gyration  in  the  M-coordinate. 

(2)  Usually,  the  momentum  (moment  of  momentum) 
equation  for  predicting  the  cylinder’s  translation  velocity 
(orientation)  is  represented  in  the  E-coordinate  (M- 
coordinate)  system.  Transformations  among  the  three 
coordinate  systems  are  used  to  convert  the  forcing  terms 
into  E-coordinate  (M-coordinate)  for  the  momentum 
(moment  of  momentum)  equation.  A  recursive  model  is 
developed  on  the  base  of  the  triple  coordinate  transform 
to  predict  the  cylinder’s  translation  and  orientation. 

(3)  The  Cylinder  Drop  Experiment  (CYDEX)  was 
conducted  at  the  Naval  Postgraduate  School  (NPS)  in 
July  2001  to  evaluate  the  recursive  model.  It  consisted  of 
dropping  cylinders  into  the  water  and  recording  the 
position  as  a  function  of  time  using  two  digital  cameras 
at  (30  Hz)  as  the  cylinders  fell  2.4  meters  to  the  pool 
bottom.  After  analyzing  the  CODEX  experimental  data, 
seven  general  trajectory  patterns  are  identified:  straight, 
slant,  spiral,  flip,  flat,  see-saw,  and  combination. 
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(4)  Model-experiment  comparison  shows  that  the  model 
well  predicts  the  cylinder’s  translation  and  orientation. 
However,  the  performance  of  the  recursive  model  for 
X  =  0  is  not  as  good  as  for  x  ^  0  . 
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